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Abstract 

We review some of the basic features of the Kazakov-Migdal model of induced 
QCD. We emphasize the role of Zfq symmetry in determining the observable prop- 
erties of the model and also argue that it can be broken explicitly without ruining 
the solvability of induced QCD in the infinite N limit. We outline the sort of critical 
behavior which the master field must have in order that the model is still solvable. 
We also review some aspects of the D = 1 version of the model where the partition 
function can be obtained analytically. 



I. INDUCING QCD 

Quantum chromodynamics is presently accepted as the only viable theory of the strong 
interactions. It describes many of the quantitative features of the interactions of hadrons 
at high energies. However, at energies lower than the hadronic scale it gives only qualita- 
tive information. Part of the reason is that perturbative QCD has only one dimensionless 
constant, the gauge coupling which is changed by renormalization from a dimenisonless con- 
stant into a mass scale. The perturbative regime is where all external momenta of Feynman 
diagrams are greater than this mass scale. It is thus impossible to address the infrared struc- 
ture of QCD, such as details of its spectrum and low energy interactions using conventional 
perturbation theory. 
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Another parameter of QCD which could be varied is the number of colors of quarks. 
QCD is known to simplify somewhat in the limit where the number of colors, N, is large 
0. In this limit, only planar graphs contribute to scattering amplitudes which consequently 
exhibit some of the qualitative features of the strong interactions. However, so far no explicit 
solution of QCD in dimensions greater than 2 is available in the large N limit and it has 
thus led to very few quantitative results. 

There are some features of the large N limit which are particularly appealing. First, there 
is every indication that the large N limit is quite smooth and that it exhibits confinement 
and dynamical chiral symmetry breaking which is found in the actual QCD for N=3. Also, 
1/N 2 = 1/9 is an expansion parameter which is less than one, so some large N results should 
have reasonable accuracy. Furthermore, the parameter 1/N is not dependent on the mass 
scale, so once a solution is found it is good at both high and low energies. This might also 
shed some light on the relationship between lattice and continuum QCD. 

Recently, Kazakov and Migdal have proposed a novel approach to the large N limit 
which has the hope of providing an exact solution of infinite iV QCD. They consider induced 
QCD which is obtained by integrating over the scalar fields in the lattice gauge theory with 
the partition function 



'KM 



fd(j>[dU}exp j-A^TrV[0(x)] + N J2 T±<!>(.x)U(xy)<j>(y)rf(xy)\ (1) 



<x,y> I 



where <f>(x) are NxN Hermitean matrices which reside on lattice sites x, U(xy) are unitary 
NxN matrices which reside on links < xy > between neighboring sites x and y, d<p is the 
Euclidean integration measure for Hermitean matrices, [dU] is the invariant Haar measure 
for integration over the unitary group U(N) and V[(f>] is a potential for the scalars. This 
model is invariant under the gauge transformations 

4>(x) — > uj{x)(j){x)uo\x) (2) 
U (xy) -► u(x)U(xy)u\y) (3) 

where u(x) is an element of U(N). By restricting the trace of <p to zero and the determinant 
of U to one in ([I]) we could also consider a model with SU(N) gauge symmetry. 

The partition function of the Kazakov- Migdal model ([I]) can be regarded as the 1/g 2 —>■ 
oo limit of lattice scalar QCD with action 

Z = J d<f>[dU] exp -N T±V[<f>(x)] + N £ Tr<f>(x)U(xy)<f>(y)U*(xy) + 



<x,y> 



+^E(w( n ) + ™ f ( n ))) 

9 □ / 



(4) 



This action differs from that of the Kazakov-Migdal model by the addition of the Wilson 
term, 

^(TrC/(n) +TrC/ t (D)) (5) 
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The Wilson term is the trace of a product of link operators around an elementary plaquette 
of the lattice. This term is the naive latticization of the continuum Yang-Mills action, 
TrF^ u F^ u . In (f|) □ denotes plaquettes of the lattice and U(0) a product of the link operators 
on the links on the boundary of □. 

Asymptotic freedom implies that the continuum limit of the lattice theory (Q) is otained 
by taking the bare coupling constant to zero, 

J2 - 00 (6) 

In fact, if instead of a lattice cutoff we had a large momentum cutoff A the bare coupling 
which would be necessary to insure one loop renormalizability of QCD is 

p^HA'/tf (7) 

The hypothesis of the Kazakov-Migdal model is that the scalar QCD might still find a way 
to arrange things so that a continuum limit exists in the opposite limit, where the bare 
coupling constant is infinite 

^0 (8, 

in || they present a naive argument to show how this might be possible. They begin with 
QCD coupled to scalars and without a kinetic term for the gluon field. The Yang- Mills 
action is induced by the vacuum polarization of the scalar fields in the cutoff theory. The 
one-loop result is 

ln(A 2 /m 2 ) (9) 




9Ld 967r : 

where m is the scalar mass. This can produce the g 2 in necessary to obtain a continuum 
limit for the gauge field sector of the theory with the ultraviolet cutoff replaced by the scalar 
mass if we take the mass of the scalar to be 

A2 \ 1/21 

(10) 

In this way, by giving up on finiteness of the scalar mass we can, at least at one-loop order, 
induce a renormalizable action for QCD. Of course this is only a rough argument. Higher 
order corrections from hard gluons will change this result significantly. They can only be 
compensated by some strong self-interactions of the scalar field. The resulting picture is 
one of a complicated, strongly interacting theory. It also requires that we have the ability 
to arrange that the scalar mass goes to infinity with a slower exponent than the cutoff in 
the continuum limit. This is possible if the scalar field theory has a socond order phase 
transition and the accompanying critical behavior. The appealing feature of this model is 
that one may be able to solve it in the large N limit. 

Recently we have proposed a slight modification of this idea |]|. The solvability of the 
large N limit comes about through the absence of a kinetic term for the gauge fields in the 
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bare Lagrangian. To leading order in N this property is also there if the kinetic term is not 
zero but is sub-leading in large N, or 

7-Jt 

where A ~ 1 and 1/g 2 ~ 1/N as N — > oo In [p]] we argued that, by tuning A appropriately 
we could still produce QCD with a string tension which is finite in the continuum limit. 

The key to the solvability of the Kazakov-Migdal model is the fact that the single-link 
Itzykson-Zuber integral can be done analytically |25l , p6| 



iiz 



[dU ]e N ^M = (12) 



where 



A[0] = det(0 i ) i - 1 = II(^-^) ( 13 ) 

13 i<j 

is the Vandermonde determinant for <ft. 

This allows us to express the partition function ([!]) as an integral over the eigenvalues of 
the scalar field 6 at each site 
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x ,i <xy> A[(f)(x)]A[(J)(y)\ 



Here, the integral over U matrices can be obtained explicitly only because the Wilson 
kinetic term is absent from the action. If the action has a Wilson term with coefficient A 
which is of order one in the infinite N limit, the Wilson term is of order TV whereas all other 
terms in the action are of order N 2 . Then, the Wilson term can be ignored and the effective 
action for the scalar field, to leading order in N is still given by flip). 

The eigenvalues <pi behave like a master field since the large N limit in fll4|) is the 
classical limit and the integral can be performed by saddle point approximation. Migdal Q 
has derived integral equations which are obeyed by the eigenvalue density and has given an 
expression for the asymptotics of the solution. Corrections to the classical behavior and the 
spectrum of elementary excitations can also be computed || . (This is of course neglecting 
the corrections which would arise from the presence of the Wilson term, which should begin 



to contribute at this order.) This model has been considered further in - |22 |. 

If the model (ffl) has a second order phase transition and if the fluctuations in the vicinity 
of the critical point are non-Gaussian, one might expect that the critical behavior should 
be represented by QCD, the only known nontrivial four dimensional field theory with non- 
Abelian gauge symmetry. Since the rough argument leading to (flO|) indicates that the 
scalar mass should scale to infinity slower than the lattice scale in the continuum limit, it 
is necessary that the scalar field exhibits critical behavior, i.e. that the effective scalar field 
theory in (|14"D has a second order phase transition itself. 

In order to familiarize the reader with the Kazakov-Migdal model, we shall begin by 
reviewing how it can be solved in the simple case of a lattice with a single site and subse- 
quently the case where the lattice is one-dimensional and periodic. These simple models are 
interesting in that one can obtain critical behavior for the scalar fields when the scalar field 
action is quadratic. 
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II. KAZAKOV-MIGDAL MODEL ON A SINGLE SITE 



It is instructive to consider the Kazakov-Migdal model on a single site. To this end 
consider the following integral 

Z = jv<j)VU e -^Tr(^)+Tr (WC /-i) (15) 

The evaluation of this integral follows very closely the method of D'Adda et. al. [ID|. The 
integral over U can be done using the Itzykson-Zuber formula |?5] and the result 
depends only on the eigenvalues of 0. We thus have that 

Z oc J A 2 (0) e-^E.*? A V 2(0) J (16) 

The determinant can be written explicitly as a sum over permutations 

Z oc £ e(o-) /n#i e- m2 Ei^ e E 4 ^W (17) 

where e(cr) is the sign of the permutation a. This integral is a Gaussian integral which can 
be done explicitly. The integrand can be simplified by introducing two real quantities a and 
b such that 



so that 



and 



a 2 + b 2 = m 2 and ab = 1 (18) 



Z oc £ e(<r) /ij^i e-K^+^w) 3 (19) 



0,6 = m 2 ± Vm 4 - 1 (20) 



Note that Eq. ( |19"D is symmetric under the interchange of a and b. 

The calculation proceeds by a change of variables from the eigenvalues fa to 

& = acj)i + 60 CT(i) (21) 

The transformation is linear and the Jacobian in independent of <fi. To evaluate the Jacobian 
first note that any permutation a can be written as a product of r-cycles of the form 

c 1 ,c 2 ...c r -> c 2 , c 3 ...c r , ci (22) 

Consider a fixed permutation cr and write a as a product of rti cycles of size 1, n 2 cycles of 
size 2, ... nk cycles of size k, ... . Clearly 

N 

£^ = N (23) 
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It is also straightforward to prove by recursion that the Jacobian for an r-cycles is 

\ 



det 



/ a -b 
a -b 0... 



... a -b 
\-b ... a ) 



a r - b r 



(24) 



where the above is an r x r matrix. 

Using the change of variables in Eq. ( f0]| ) we find that 

z oc e ( a ) Z 

with 



(25) 



n« 



<V 2 



n n 

cycles i£ cycles 



1 



a r - b r 



(26) 



Performing the Gaussian integral we see that 



z a = Tr^n 



N 



i=l 



a* - 6 l 



oo 

c* II 

i=l 



a i - 6* 



(27) 



since nfc=0 for k>N. 

The next step is to do the sum over permutations fl25|). To this end note that there are 



N 



N \ n 



1 1 



(28) 



distinct permutations which can be described the a given set (n 1; ...n N , ...) (with nfc=0 for 
k>N). Thus 



OO CO OO OO / 1 

z oc £ ••• E •••*(Jv-E^fc) n 

m=0 ^fc=0 fc=i 
where we have used the fact that 



_L JL (29) 



(30) 



k=l 



Using the integral representation of the 5 function 5(x) = ^ d9exp(— ix9) we find that 



oo / oo AQkn 



z oc de e- im n E 



£t x [k(a k - b k )] nk n k \ 



(31) 



The sum in Eq. ( |31| ) can be summed to an exponential. The remaining product then 
becomes a sum over the exponents so that 
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Z cc I dO e 



-iNO 



cxp 



E 

k=l 



'-ae i6 ) k 



1 



,2k 



(32) 



k 1 

where we use have used the fact that b—l/a. Without loss of generality we assume that a<l 
(i.e. a=m 2 — a/ m 4 — 1 . The term 1/(1 — a 2k ) can now be expanded in a geometric series 
and the resulting sum over k sums to a logarithm. Thus 



Z oc / d0 e - *™ J] exp [log (l + e l 

n=0 



i6 n 2n+l 



dd e 



-iNO 



n 

n=0 



a 2n+1 



(33) 



The final step is to perform the integration over 9. To do this note that only terms in 
the product which are proportional to exp(iN6) will give a nonzero integral. Thus 



ni=0 7i2=ni+l nN— n N+i 



(34) 



The sums can be done one at a time. They are all geometric series. The result is 



JV 



,2\ k 



Z « TvII 



il 



/ 2 N^/2 " 1 



(35) 



In fact if we define q=a 2 = (m 2 — -y/m 4 — l) 2 then 



Z (x q 



N 2 /2 



N 



Vfc=l 



— 



(36) 



Notice that apart from the factor g^ 2 / 2 this is just the expression for a q-factorial. It is 
interesting to note that a similar expression is obtained for the partition function of a system 
bosons on a circle at nonzero temperature. In fact the partition function for a bosonic string 
on a circle (suitably restricted to the singlet sector) is precisely given by the above formula. 



III. KAZAKOV-MIGDAL MODEL ON A CIRCLE 

It is straightforward to generalize the above calculation to the evaluation of the partition 
function of the Kazakov-Migdal Model on a circle. This calculation is discussed in Ref. [T^ 
and 0. We consider the partition function 

]J T><f> x n T>U X>X+1 e~ m2 Tr Tr (*- u **+i**+i u £ + i) (37) 

2=1 X=l 

where the fields (p x live on the sites x of a circle with periodic boundary conditions (pL+i = 4>i 
and the U XjX+ \ live on the links of the circle. There are two distinct ways of calculating this 
partition function. The first method is to eliminate almost all the U x>x+ iS by a gauge 
transformation leaving only one U which cannot be eliminated. This can be chosen to 
be the U on the first link i.e. U=Ui^- The result is a partition function which involves 
an integral over only a single link variable, and a Gaussian integral over all N Hermetian 
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matrices. The second method is to integrate over all the link variables U XiX +i explicitly 
using the same formula as for the single link case above. This second method generalizes 
the method presented in the previous section for the single link integral. The details of this 
calculation are discussed by D'Adda and Panzeri | 10[| . The basic idea is to perform all the 
U integrations in Eq. fl3"T|). The result is 

r L ( N \ ( L det e^^+i \ 

z "/l(S*) 4!w (S») (38) 

where 4> l x is the i'th eigenvalue of the matrix <f> x . Note that the Vandermonde determinants 
A(0) precisely cancel leaving us with 



Z oc 



/ n (^t e+W*» e~ m2 ^M) 2 \ (39) 

x=l 



This is now a Gaussian integral of precisely the same form as the integral for the K-M model 
on a point. It is evaluated in Ref. fl(|. The result is 



, N 1 

ft In I TT 1 



Kk=l * 'I 

with 



z«9 LN > 2 [Urzhs («) 



q = m 2 - y/m 4 - 1 (41) 



Note that this is obtained from the single link integral by a simple replacement of q by q 
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IV. SYMMETRY AND OBSERVABLES 

It was pointed out in that, like all adjoint lattice models, the Kazakov-Migdal model 
(H) has an extra gauge symmetry which is not a symmetry of continuum QCD. The action 
in ([!]) is invariant under redefining any of the gauge matrices by an element of the center of 
the gauge group, 

U(xy) — ► z(xy)U(xy) (42) 
<j){x) -> <j>{x) (43) 

where z{xy) G U(l) if the gauge group is U(N) and z(xy) G Z N if the gauge group is SU(N). 
(We shall call the symmetry a Zn gauge symmetry in either case.) It was subsequently 



pointed out by Gross [11| and by Boulatov [19| that there is a larger symmetry of this kind: 
one could redefine U and W by any element which commutes with the matrix (p. The Zn 
symmetry in ( ^3|) is the maximal subgroup of the transformations discussed by Gross and 
Boulatov which can be implemented with field independent elements z(xy). Because of this 
symmetry the conventional Wilson loop observables of lattice gauge theory have vanishing 
average unless they have either equal numbers of U and operators on each link or else, 
in the case of SU(N), unless they have an integer multiple of N U's or iV W's. 
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In conventional QCD, the expectation value of the Wilson loop operator gives the free 
energy for a process which creates a heavy quark-antiquark pair, separates them for some 
time and lets them annihilate. From the asymptotics for large loops, one extracts the 
interaction potential for the quarks. If the expectation value of the Wilson loop behaves 
asymptotically like e~ aA where A is the area of a minimal surface whose boundary is the loop, 
the quark-antiquark potential grows linearly with separation at large distances and quarks 
are confined. The parameter a is the string tension. On the other hand if the expectation 
value of the Wilson loop goes like the exponential its perimeter then the potential is not 
confining. 

In the Kazakov-Migdal model (JJ), due to the Z N symmetry, the expectation value of the 
Wilson Loop is identically zero for all loops with non-zero area. We can interpret this as 
giving an area law with infinite string tension, a = oo, and no propagation of colored objects 
is allowed at all. (An exception is the baryon (U N ) loops in the case of SU(N) where the 
correct statement is that iV-ality cannot propagate.) It is for this reason that the original 
Kazakov-Migdal model has difficulty describing pure gluo-dynamics. 

The Zn symmetry of the pure Kazakov-Migdal model is broken explicitly by the intro- 
duction of a Wilson term in (|4]). However, if the Wilson term has vanishingly small coefficient 
and is negligible in the large N limit, one might expect the problem of Z^ symmetry to 
remain - the string tension would still be infinite. It has been argued in that this need 
not be so. It was shown that, if the scalar fields exhibit a particular kind of critical behavior, 
it is possible that the presence of an infinitesimally small Wilson term is sufficient to give a 
finite string tension. 

There are currently several other points of view on how to avoid the constraints of 
Z N symmetry. In |J it was suggested that if there is a phase transition so that the Z N 
symmetry is represented in a Higgs phase, the resulting large distance theory would resemble 
conventional QCD. This approach has been pursued in ||, [12|, ||13|| . An alternative, which 
was advocated in [§,0, is to use unconventional observables such as filled Wilson loops which 
reduce to the usual Wilson loop in the naive continuum limit but which are invariant under 
Z N . The third possibility is to break the Z N symmetry explicitly. This was suggested by 
Migdal |15j in his mixed model in which he breaks the Z^ symmetry by introducing into 
the model heavy quarks in the fundamental representation of the gauge group. 

In []]J we considered explicit Zn symmetry breaking using a Wilson term. We showed 
that it is equivalent to using the filled Wilson loop observables. We used that fact that the 
filled Wilson loops arise naturally from ordinary Wilson loops in a modified version of the 
Kazakov-Migdal model which has additional explicit symmetry breaking terms. We argued 
that one version of this modified model should be solvable in the large N limit. 



V. FILLED WILSON LOOPS 

We begin with a brief review of the properties of the filled Wilson loop operators which 
were introduced in |J and discussed in detail in |7[]. These are a special class of correlation 
functions which survive the Zn symmetry of the original Kazakov-Migdal model ([I]). They 
are defined by considering an oriented closed curve T made of links of the lattice. The 
ordinary Wilson loop operator on T is given by 
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W[T]=Trl [] U(xy)\. (44) 
[<xy>er J 

For any surface S which is made of plaquettes such that the boundary of S is the curve T 
we define 

W F [T, S] = W[T] J] W*[n] (45) 
aes 

where □ denotes an elementary plaquette in the surface S. The filled Wilson loop for T is 
now defined as 

W F [T}=J2KS)W F [T,S] (46) 

s 

where the sum is over all surfaces S whose boundary is the loop T with some (yet to be 
specified) weight function /i(S'). Notice that for each plaquette □ G S we have inserted the 
negatively oriented Wilson loop Thus for arbitrary weight functional n{S) the filled 

Wilson loop operator is invariant under the local gauge symmetry since it has equal 
numbers of U and operators on each link. Although we have assumed that the loop is 
filled with elementary plaquettes this can be easily generalized to other fillings (the other 
extreme case being the adjoint loop [r]P0[T]). We can also define the 'filled correlator" 
of more than one loop by summing over all surfaces whose boundary is given by those loops. 

In ref. it was shown that computing the expectation value of W^fT] is equivalent to 
computing the partition function of a certain statistical model on a random two-dimensional 
lattice. When computing Zn gauge invariant correlation functions of [/-matrices in the 
master field approximation the 0-integral is evaluated by substituting the master field = 
diag(0i, . . . , 07v ) for the eigenvalues of 0. 

<T r rrt , [ #[g^g ^U hjl ■ ■ ■ Uj h ■ ■ ■ 

Hn ' ' ' klh " Jd<f>[dU)e-^vm-j:w<t>ui) 



Jd4>[dU\e^^m Uixh ...U, 



t 



Jd0[dU}e Tl ^^ ut ) 

In this integral, the scalar field is written at = V^bV^ with (p£> a diagonal matrix. The 
eigenvalues of are fixed at the value of the master field and, in order to obtain gauge 
invariance of the correlator, the angular matrices V are still integrated, 

rf0 = A(4>)DV (48) 

In any gauge invariant correlator, the matrices V can be absorbed by redefining U. 

If we consider for the moment surfaces which are not self-intersecting so that the filled 
Wilson loop correlator has at most one UU* pair on any link we need to consider only the 
two field correlator < U^U^ >. Gauge invariance implies that 

<UuUl>=C ij 6 u 6 jk with Cu = J m J ^ ,_ (49) 

3 3 3 J d4>[dU]e N ^(* u * W ) 
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Thus, in the master field approximation, the expectation value of the filled Wilson loop is 
given by 



<w f [t] >=EMS) 

s 



n e n c n*) 

. sites i(x) = l J links 

\xeS I <xy>es 



(50) 



This is a generalized Potts model on a random surface in which iV-component spins reside 
at each site and the Boltzmann weights for the bonds are correctly normalized to be 
conditional probabilities; Y,iCij=l, J2j Cij=l. 

VI. COMPUTING CORRELATION FUNCTIONS 



Techniques for evaluating CL- for general N and for arbitrary <p are presented in [17| and 



|Tj|. An explicit formula for SU{2) is given in Although the general formula for Cy in 
SU (N) is quite difficult to deal with, it is still possible to estimate the surface dependence 
of the statistical model partition function in ( |5D| ) when the master field is homogeneous 
by considering two different limits. 



A. Low Temperature 

First, consider the case where 0« are large. We also assume that the eigenvalues <pi are 
not too close to each other in the sense that 

E T- l -N2 « iV ( 51 ) 

#j [4>i - <Pj 

(Note that this is a single sum over j for fixed i.). The integral in ([12]) is known to be 
exact in the semi-classical approximation (see |7]] for a discussion). The classical equation 
of motion is 

[U(j)rf,<t>] = (52) 

which, since is diagonal, is solved by any U of the form Uq = DP where D is a diagonal 
unitary matrix and P is a matrix which permutes the eigenvalues, 

(P$P)ij = (%0p W (53) 

Also, when N is large and is not too small the identity permutation gives the smallest 
contribution to the action in (12) and therefore is the dominant classical solution. In this 



case we use this minimum to evaluate the correlators, 

Jfi / [dU}e N ^^U hjl . . . U injn Ul lh ■ ■ ■ Uli n = <U ■ ■ ■ *Wn5&& (54) 
where we have written the normalized integral over diagonal matrices 
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si\;x = iUede e u P<q ^(6 p - ^)e i (^ + - + ^-^--^)//n £ ^n P<g sin 2 (^ - 6 q ) 

1 if ii . . . i n is a permutation of k\ . . . k n , , 

otherwise 

We have decomposed the integration over unitary matrices into an integration over the 
diagonals and an integration over the unitary group modulo diagonals |27[| . The diagonals are 
the 'zero modes' for the semiclassical integral and must be integrated exactly. The unitary 
modulo diagonal integral is damped by the integrand and is performed by substituting the 
classical configuration. Of course, to get the next to leading order the latter integration 
must be done in a Gaussian approximation. It can be done for the first few correlators, The 
result is '| '| 

= 1 — T7 E — 2 w hen i ^ j Cij = -r;— — f 

N k+i ((pi - ■ N ^. _ 0^ 

We remark that similar calculations can be easily done for correlators of more than two 

C/'s 

CaM = !Tz J [dU^^^lU^lU^ 2 (56) 



Cu t kk — i — j7 E 



+ 



i*5 , k = i CiS ^NW^W (57) 

where the next corrections will be of the form 

1 1 



E -r — -T2 ■ (58) 



N2 ((pi - m+k 

We call this limit of large the "Low Temperature" limit since in this limit 

C£ = 5 lj + ... (59) 

and the value of the spin at each site is equal. In this case the represent the Bolzman 
weight for a perfectly ordered system. These two cases lead to profoundly different behavior 
for the filled Wilson loop. We shall assume that, by choosing the potential for the scalar 
field in (|14]) appropriately, either of these limits could be obtained (the eigenvalue repulsion 
due to the Vandermonde determinants in fll4]) and the possibility of adding repulsive central 
potentials makes the low temperature limit more natural). 
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B. High Temperature 



The other limit is where <\> is small. There, we can obtain the correlators by Taylor 
expansion, 

C l3 = J [dU] (l + NTt$U4rf + ■ ■ •) \Uij\ 2 = ^ + ^ + • • • (60) 

We call this the "High Temperature" limit since in this limit is independent of i and j. 
It thus represents the Bolzman weights for a highly disordered system. 



C. Renormalization of the String Tension 

We begin by estimating the value of the filled Wilson loop for a fixed surface S. In the 
"high temperature" case the statistical model is disordered. The sums over configurations 
at the various sites are independent and they contribute an overall factor N v (where V is 
the number of vertices on the surface) to the expectation value of the filled Wilson loop. 
Furthermore each link contributes a factor Cij—l/N so that the links contribute a total 
factor of N~ L where L is the total number of links. It follows that the expectation value of 
the filled Wilson loop goes like 

< W F [T, S] > HT ~ N y - L = N 2-MS) N -A { s) (61) 

where A(S) is the area and g(S) is the genus of the surface S (i.e. the number of plaquettes 
comprising S) and we have used Euler's theorem, % = 2 — 2g = V — L + A. We thus get the 
renormalization of the string tension 5«ht — log N. Notice also that higher genus surfaces 
are suppressed and that the loop (genus) expansion parameter is 1/N 2 . This is precisely 
what is obtained in the conventional strong coupling expansion of Wilson's lattice gauge 
theory which is known to describe a string theory with extra degrees of freedom associated 



with self-intersections of the string [g3 



In the "low temperature" case, the statistical system is ordered. The spins on all the 
sites are frozen at a uniform value. In this case the partition function is proportional to the 
degeneracy of the ground state, 

<W F [T,S}> LT = N (62) 

Note that in this case the statistical model gives no contribution to the string tension 
(5ct LT « 0) and there is no suppression of higher genus surfaces. 

In order to proceed to the evaluation of the filled Wilson loop we need to choose a weight 
function fi(S) in order to perform the sum over surfaces. The most reasonable criterion for 
choosing such a weight function is our desire to get a finite physical string tension in the 
continuum limit. In order to accomplish this goal we must choose a weight function fJ^(S) 
which depend on the area of the surface differently in the low and in the high temperature 
cases. It is known that the number of closed surfaces with a given area grows exponentially 
as 

n{A) ~ A K{9) e^ A (63) 
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where n{g) is a universal constant which depends only on the genus of the surface and /xo is a 
non- universal, regulator dependent constant which will lead to a renormalization of the 
string tension. In our case, although the surfaces are open, the above formula should still be 
valid for surfaces whose area is much larger than the area of the minimal surface bounded 
by T. If the continuum limit of our theory is realized in the "high temperature" phase we 
should use the weight function (jl ht (S) ~ ]\[ A (s) e -^oA(S) _ This leads to a vanishing string 
tension in the lattice theory which is a necessary condition for having a finite string tension 
in the continuum limit. To accomplish the same goal in the "low temperature" phase we 
should use fi LT (S) ~ e -voA(S)_ Although these choices of fi(S) give the desired result, it is 
rather unnatural to have to choose in such an ad hoc fashion. 



D. Filled Wilson Loops from the Wilson Action 



Fortunately there is a very natural way to obtain the sum over surfaces in fl5U|). Consider 
the following expectation value 

< W F [T] >= [ ft 7 , (64) 

where W[T] is the conventional Wilson loop. Remember that the average is weighted by the 
Kazakov-Migdal action as in ([47]): In the master field limit it is computed by integrating only 
over [/-matrices with = and with the Kazakov-Migdal action. Note that the exponent 
in (B3) is simply the conventional Wilson kinetic term for the gauge fields in lattice gauge 



theory. If we expand the right hand side of (p4f) in A the non-vanishing terms are all of those 



surfaces which fill the Wilson loop. The result is thus a filled Wilson loop with a surface 
weight fi(S) = X A{S) . 

It is clear that we would obtain exactly the same expression (in the master field ap- 
proximation) by evaluating the expectation value of the ordinary Wilson loop operator in 
the modified version of the Kazakov-Migdal model (^) in which a conventional Wilson term 
(A 2d + W^p])) is added to the action. This term breaks the Z N gauge symmetry 

explicitly and allows Wilson loop operators with non-zero area to have non-zero expecta- 
tion values. We would expect that it is necessary to keep A small if one is to maintain the 
successes of the Kazakov-Migdal model. We shall now argue that in the "low temperature" 
limit this picture is self-consistent in the sense that the physical string tension is finite when 
A is small and consequently the saddle point solution of the original model is unchanged. 
We shall also see that this is not the case in the "high temperature" phase. 

Let us begin by determining how A should behave in the continuum limit if we are to 
have a finite physical string tension. As discussed above a necessary condition for having 
a finite physical string tension is that the string tension in lattice units should vanish. It 
is thus necessary for the bare string tension — In A to be chosen so as to precisely cancel 
the renormalization of the string tension due to both the statistical model to the sum over 
surfaces. It is straightforward to check that in the "high temperature" phase we must choose 
Aht = iVe _M0 , whereas in the "low temperature" phase we must chooseA LT = e~ M0 . Notice 
that this A is proportional to N in the "high temperature" phase and thus cannot be assumed 
small in large N. 
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In the large N limit of conventional lattice gauge theory the coefficient of the Wilson 
term must be proportional to N if one is to obtain a consistent large N expansion. In our 
case we see that this is true for the "high temperature" phase in which case the Wilson term 
is of the same order as the Kazakov-Migdal term and it thus plays an important role in the 
infinite N limit. One can say that in this phase we have ordinary QCD. Unfortunately it is 
impossible to preserve the master field solution of the Kazakov-Migdal model in this limit 
since the Wilson term, being of order N, would modify the large N solution, ruining the 
self-consistency of the mean-field approximation as described here. 

The situation is much more appealing in the "low temperature" phase. In this case the 
required coefficient of the Wilson term is of order one. It is subdominant and therefore 
negligible in the large N limit. Thus, Migdal's solution [|J of the Kazakov-Migdal model 
in the large N limit should still apply to our proposed modification of the action. In fact 
the only reason that the Wilson term is important at all in the large N limit of the "low 
temperature" phase is related to the collective phenomenon which orders the statistical 
system on the surfaces. It effectively makes the statistical model's contribution to the string 
tension much smaller than would be expected from naive counting of powers of 1/N and a 
truly infinitesimal breaking of the Zn gauge symmetry (A LT /iV — > as N —>■ oo) is sufficient 
to make the averages of Wilson loop operators non-vanishing. The self-consistency of this 
picture can also be demonstrated by computing the contribution of the Wilson term to the 
free energy. This can be computed in a small A expansion. For a cubic lattice the result is: 

Z =< ^(EnW^M)) >= Zkm e*JNV D{D ~ 1} (AL + 2A 6 LT + . . )) (65) 

is of order N (where V is the volume, D is the dimension). This should be compared with 
the free energy in the pure Kazakov-Migdal model which is proportional to N 2 . Here, the 
first term in the free energy is the contribution of the doubled elementary plaquette and 
the second term is due to the two orientations of the elementary cube. It is interesting 
that, to order 6, there is no energy of interaction of doubled elementary plaquettes with 
each other. We conjecture that the interaction energy of surfaces is absent to all orders and 
the free energy obtains contributions from all possible topologically distinct surfaces which 
can be built from elementary plaquettes. This suggests a free string picture of the "low 
temperature" limit of the Kazakov-Migdal model at lattice scales. 

E. Self-Intersecting Surfaces 

We have thus far neglected the self-intersecting surfaces in the sum (|50| ) which are gen- 
erated by the expansion of ( |6~4l) in A. In order to evaluate the contribution of these surfaces 
we need to compute the correlator of n UU^ pairs on the same link. The computation of 
these correlators in full generality is quite complicated. In the Appendix we compute them 
in the "low temperature" (ordered) phase. We find that 

< U hh . . . U lnJn Ut lh ■ ■ ■ Ul nln >= 5 lin . . . 5 knln S£± (66) 

where S^'" % j^ is the tensor which is one if i\ . . . i n is a permutation of k± . . . k n and is zero 
otherwise. It is now evident that in this limit the [/-matrices are replaced by unit matrices 
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which freeze together the spin degrees of freedom on the various intersecting surfaces. As a 
special case we can consider a single, connected, self-intersecting surface. In this case all the 
spin indices on the surface are equal and since = 1 when all arguments are equal the 
partition function of the statistical model corresponding to that surface is simply N just as 
it was for a non-intersecting surface. Thus just as the statistical model does not contribute 
to the string tension it also does not contribute to the interaction energy of self-intersecting 
surfaces. This implies that in the "low temperature" limit, the sum over connected surfaces 
which have a common boundary behaves like a Nambu-Goto string theory with no internal 
degrees of freedom. 



VII. DISCUSSION 

In summary, the self-consistency of the "low temperature" limit leads us to a new large 
N limit of the conventional lattice gauge theory coupled to scalars: 

d(p[dU] exp I —N^2TiV[(p(x)] + N £ Tr<f){x)U{xy)<f){y)U*{xy)+ 

\ x <x,y> 

+A^(D) + ^(0))) (67) 

The conventional large N limit occurs when A is of order N and describes scalar QCD. 
The other limit occurs when N — > oo with A of order one. This model is soluble using the 
Kazakov-Migdal approach. 

It is the latter case in which A remains constant that is of special interest to us. In this 
case we saw that the large N expansion corresponds to a string theory with some unusual 
features. The partition function and the Wilson loop expectation value can be described as 
a sum over surfaces. What is unusual is that the genus of the surfaces is not suppressed in 
the large N limit, as it is in continuum QCD. (We do of course expect the higher corrections 
in 1/N to suppress higher genus terms.) For a continuum string theory this sum over the 
genus is badly divergent. This, together with the presence of tachyons, suggests that the true 
ground state of the string theory is some sort of condensate. This could pose a complication 
for the present version of the Kazakov-Migdal model in the continuum limit and deserves 
further attention. It is still a mystery to us how the sum over all surfaces at the lattice scale 
should turn into the sum over planar diagrams in the continuum theory of QCD. 

An alternative to the model presented here is the mixed model which was invented 



by Migdal [I5j to solve the problem of Zn symmetry. It contains heavy quarks in the 
fundamental representation of the gauge group. Despite the obvious differences between our 
model and Migdal's mixed model they have many features in common. As in all cases when 
there are fields in the fundamental representation, the asymptotics of the Wilson loops in 
the mixed model exhibit a perimeter law. In conventional QCD one would expect that if 
the quarks are heavy enough, there is an area law for small enough loops, i.e. there would 
exist a size scale which is far enough into the infrared region that the quark potential is 
linear but the interaction energy is not yet large enough that it is screened by producing 
quark- ant iquark pairs. Thus, in QCD we expect that adding heavy quarks would not ruin 
the area law for Wilson loops smaller then some scale. 
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The mixed model has just the opposite scenario, it is possible to get an area law only when 
the heavy quarks are light enough. This is a result of the fact that, in the Kazakov-Migdal 
model, no Wilson loops are allowed at all unless the Z N symmetry is explicitly broken. In 
the mixed model, the Z N charge of links in a Wilson loop must be screened by the heavy 
quarks. This can happen in two ways. First, the Wilson loop can just bind a heavy quark 
to form an adjoint loop - giving a perimeter law for the free energy of the loop. This is 
the leading behavior if the fermion mass, M, is large. The free energy would go like 1/M P 
where P is the perimeter. The only way an area law might arise is when the fermions are 
light enough that their propagators could from a filled Wilson loop with free energy 1/M 2L 
where L ps 2A is the number of links. Then, since the entropy for filled loops is much larger 
than that for adjoint loops, these configurations would be important if M 4 < e M0 . Then, the 
asymptotics behavior of the Wilson loop would still have a perimeter law but there would 
be loops with 4A — P < // /lnM where there would be an approximate area law. 
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